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Abstract 

The inclusive charm hadron semileptonic branching fractions B(c — > e) and B(c — > //) in 
— > cc events have been determined using 4.4 million hadronic decays collected with 
the OPAL detector at LEP. A charm-enriched sample is obtained by selecting events with 
reconstructed D*^ mesons. Using leptons found in the hemisphere opposite that of the 
D*^ mesons, the semileptonic branching fractions of charm hadrons are measured to be 

B(c ^ e) = 0.103 ±0.009 and B(c ^ /i) = 0.090 ± 0.007 , 

where the first errors are in each case statistical and the second systematic. Combining 
these measurements, an inclusive semileptonic branching fraction of charm hadrons of 



B(c -^i) = 0.095 ± 0.006 t^oml 



is obtained. 
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1 Introduction 



The inclusive charm hadron semileptonic branching fractions B(c — e) and B(c — » /i) are 
defined as the average of the semileptonic branching ratios of weakly decaying charm hadrons 
weighted by their production rates in prompt charm events, Z° —* cc. Inclusive semileptonic 
branching ratios are a means to investigate the dynamics of heavy quark decays, and have been 
studied in much detail for bottom quarks j^. The inclusive semileptonic branching ratio of 
charm hadrons has not previously been measured at LEP, even though it is an important input 
to a number of measurements performed at energies around the Z° resonance . 

The inclusive semileptonic branching fraction of charm hadrons has so far been measured 
at centre-of-mass energies significantly below the Z° mass 0,^. Many of these measurements 
depend strongly on the modelling of the h —>■ i background in the sample. In this paper a 
measurement of B(c — > e) and B(c — >■ /x) is presented which is much less dependent on the 
bottom background, since it is done in a sample of events enriched in Z° cc decays. This 
sample is prepared by selecting highly energetic D*+ mesonsQ. The hemisphere opposite to the 
one containing the D*+ meson is searched for a lepton, yielding a measurement of the inclusive 
semileptonic branching fraction of charm hadrons. 

The paper is organised as follows. The principle of the analysis, in particular the method 
used to subtract the background, is discussed in section ^. After a brief review of the event 
selection in section ^, the identification of charm events using reconstructed D*+ mesons is 
described and the determination of the charm fraction in the sample is summarised in section 
^ The preparation of the lepton sample in charm-tagged events and the measurement of the 
background in this sample is described in section ^, followed by the presentation of the results 
in section |^. Systematic errors are given in section |^. 



2 Analysis Principle 

A sample of Z° cc enriched events is found using reconstructed D*+ mesons. Each event 
is divided into two hemispheres by the plane perpendicular to the thrust axis. Leptons are 
searched for in the hemisphere opposite the D*"*" meson. Background is suppressed by requiring 
that the D* and the i have opposite charge. The number of leptons found in the hemisphere 
opposite that of the D* meson has contributions from prompt charm decays, c — from 
prompt bottom decays, b — >■ from cascade decays, h ^ c ^ i, and from background. It can 
be written as 

(1 - JT") [Xeff B(b ^ i) e'r' + (1 - Xeff) B(b ^ c ^ £) e'l-'^-'] } + Ar+^ .(1) 

Here A^d*+ is the number of D*"*" mesons found in the data sample, /°*^ is the fraction of 
these D*"*" mesons coming from TP — > cc events, and is the number of background events, 
where either a D*"*", a lepton or both are misidentified, but where the charge correlation is 
correct between the two hemispheres. This background will be denoted as "combinatorial 
background" . The parameter Xefr is the effective mixing parameter for the mixture of neutral B 



^Throughout this note charge conjugation is always imphed, unless explicitly stated otherwise. 



4 



mesons selected, and e^^^, e^^^ and are the efficiencies to find a lepton opposite a D*"*" 

in the channel indicated, with the correct charge correlation. To simplify this and the following 
equations, leptons produced in b — > c — ^ £ decays and b — > r ^ £ decays are included in the 
b — >■ £ decays. Since a pair of leptons, one with the correct and one with the wrong charge 
correlation is produced in b — > J/\E' — decays they are equally split between the b — >^ £ 
and the b ^ c — > £ decays. 

The goal of this analysis is the measurement of B(c It is extracted from A^£.+ the 

number of TP cc events where simultaneously a D*"*" meson in one hemisphere and a lepton 
in the opposite hemisphere is found: 

iV^.,,,_ = iVD*+ /°*" B(c ^ I) er' . (2) 

A sample of events which does not contain contributions from prompt charm decays is prepared 
by selecting events where the D*"*" and the lepton have equal charge: 

A^D*+,.+ = iVD*+(l-/e°*^) {(1 - Xcff) B(b £) e^^^ + xeff B(b ^ c ^ £) e^-^^^^} + iV++ . (3) 

Here is the number of combinatorial background events in this wrong sign sample. The 
number of leptons from charm hadron decays can be calculated by solving the two equations 
and 1^ for A'^.^. ^_ defined in equation |^. The solution can be written in terms of the difference 
of the two samples of events and two small corrections, 

Ar^.+_,_ = (ArD*+,,- - Nj,,+^,+ ) - AATb - AA^bgd • (4) 

The first correction, AA^b? can be derived directly from equation |l] and equation |^ and refiects 
the fact that mixing affects both samples differently. It is calculable from the known branching 
ratios and the mixing parameter: 

AA^b = A^D*+ (1 - fT^) (1 - 2xeff) {B(b ^ c ^ £) e'l-''-' - B(b ^ £) e,"^^} . (5) 

The second correction, AA^bgd, is the difference between the combinatorial background term 
in both samples, AA'bgd = ~ ^bgd- This number is determined using both data and 

Monte Carlo simulations, as will be discussed in section |5.1| . Finally the inclusive semileptonic 
branching ratio of charm hadrons, B(c — £), is calculated from equation ^ as 

B(c-£)=N^.,,-— -1— ^, (6) 

where the number of events with a D*"*" meson, the number of leptons found in this sample, the 
efficiencies to reconstruct the leptons in the tagged charm sample, and AA^b and AAbgd have 
to be known. Each of these inputs will be discussed in the following sections. 

3 The OPAL Detector and Event Selection 

A detailed description of the OPAL detector can be found elsewhere p. The most relevant 
parts of the detector for this analysis are the tracking chambers, the electromagnetic calorimeter, 
and the muon chambers. The central detector provides precise measurements of the momenta 
of charged particles by the curvature of their trajectories in a solenoidal magnetic field of 
0.435 T. The electromagnetic calorimeter consists of approximately 12000 lead glass blocks. 
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which completely cover the azimuthal range up to polar angles^] of | cos6\ < 0.98. Nearly the 
entire detector is surrounded with at least three layers of muon chambers, which are placed 
behind an approximately one meter thick iron magnet flux return yoke. 

Hadronic Z° decays are selected using the number of reconstructed charged tracks and the 
energy deposited in the calorimeter, as described in p. The analysis uses an initial sample of 
4.4 million hadronic decays of the Z° collected between 1990 and 1995. 

Hadronic decays of the Z° have been simulated using the JETSET 7.4 Monte Carlo model 
with parameters tuned to the data ||^. The Monte Carlo samples are about five times larger 
than the collected data sample. Heavy quark fragmentation has been implemented using the 



model of Peterson et al. with fragmentation parameters determined from LEP data pO 
The samples have been passed through a detailed simulation of the OPAL detector 0] before 
being analysed using the same programs as for the data. Jets are reconstructed in the events 
by the JADE jet finder using the EO scheme with a cut-off parameter Xmin = 49 GeV^ |]12|. 



4 Charm Tagging 

The tagging of Z° cc events is based on the reconstruction of charged D*"^ mesons in 
five different decay channels. The identification algorithm and the method to separate the 
different sources contributing to the observed D*+ signal have been presented in a previous 



OPAL paper [13|, and will only be briefly reviewed. 



The D*^ mesons are reconstructed in the following five decay channels: 

D*+ ^ D°7r+ 

U K-TT+ , "3-prong" 
' — > K^e^Uc , "electron" 
' — > K^iU+z/^ , "muon" 
U K-tt+ttO , "satellite" 

' — > K^7r+7r^7r+ , "5-prong" . 

The muon and the electron channels are collectively referred to as "semileptonic" . No attempt 
is made to reconstruct the 7r° in the satellite channel, nor the neutrino in the two semileptonic 
channels. Electrons are identified based on their specific energy loss, dE/dx, in the central 
tracking chamber and the energy deposition in the electromagnetic calorimeter. An artificial 



neural network trained on simulated events is used to perform the selection [14|. Electrons from 



photon conversions are rejected as in |TH]. Muons are selected using matching of track segments 



of the central tracking chambers and the muon chambers, as described in |T^. The purity of 
kaons is enhanced by requiring that the dE/dx measurement of the candidate is compatible 
with that expected for a kaon. If the track combination has an invariant mass M^o within 
the limits given in table 0, the combination is accepted as a D° candidate. The combinatorial 
background is reduced by a cut on the cosine of the helicity angle, cos^*, measured between the 
direction of the D° in the laboratory frame and the direction of the kaon in the rest-frame of 
the D° candidate. Background is expected to peak at ±1 in this variable, while true D° mesons 
are uniformly distributed. These D° candidates are combined with a candidate for the pion 
in the D*"^ D^vr"*" decay. Background from bottom decays and combinatorial background 



^The OPAL coordinate system is defined as a Cartesian coordinate system, with the x-axis pointing hori- 
zontally towards the centre of the LEP ring, the z-axis in the direction of the outgoing electrons, and the y-axis 
points approximately vertically upwards. The polar angle is measured with respect to the z-axis. 
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D° decay mode 


Variable 


3-prong 


semileptonic 


satellite 


5-prong 


range 


0.4-1.0 


0.4-1.0 


0.4-1.0 


0.5-1.0 


Mdo [ GeV] 


1.79-1.94 


1.20-1.80 


1.41-1.77 


1.79-1.94 


AM [ GeV] 


0.142-0.149 


0.140-0.162 


0.141-0.151 


0.142-0.149 


cos 61* xd*+ < 0.5 


-0.8-0.8 




cos6'* a;D*+ > 0.5 


-0.9-1.0 


W^if/d. ^D.+ < 0.5 


> 0.1 




Relative abundance 


0.231 


0.121 


0.355 


0.293 


Signal /background 


3.496 


3.233 


1.223 


0.879 



Table 1: List of selection cuts used in the D*"*" reconstructions. Wj^^^ is the probability that 
the measured dE/dx value is compatible with that expected for a kaon at the measured momentum. 
This cut is only applied to the kaon candidate in the D'^ decay. The background distribution thus 
obtained is normalised to the candidate AM distribution in the range 0.18 GeV < AM < 0.20 GeV 
(0.19 GeV < AM < 0.22 GeV in the semileptonic channels). In the last two lines of the table, the 
relative abundance of each channel and the signal/background ratio is given, as measured from the 
data. 



is reduced by selecting candidates with a large scaled energy, xd*+ = -E'^'t+^Z-Ebeam- The final 
selection is made on the mass difference AM = Md*+ — M^o between the D*"*" candidate and 
the corresponding D° candidate. 

If more than one D*"^ candidate is found in an event, only one candidate is accepted accord- 
ing to the following procedure. A 3-prong decay is preferred over a semileptonic one, which 
in turn is preferred over a satellite, and a 5-prong decay is selected last. If more than one 
candidate is found within the same decay channel, the one with M^o closest to its nominal 
value of 1.864 GeV [TD (1.6 GeV for the satellite) is selected. In figure 0, the mass difference 
distributions AM = Md*+ —M^o are shown for the different channels. In total 27662 candidates 
are selected in all five channels. 

The selected sample of D*"*" candidates has contributions from: D*"*" mesons produced in 
Z° — >■ cc events (signal); D*"*" mesons produced in Z° — bb events; D*"*" mesons produced 
in events where a cc pair is produced in the splitting of a gluon; combinatorial background. 
The combinatorial background in the sample of D*+ mesons is subtracted on a statistical basis 
using an independent sample of background candidate events, selected based on a hemisphere 
mixing technique first introduced in [O]. The candidate for the pion in the D*"*" — >■ D'^vr"'" 
decay is selected in the hemisphere opposite to the rest of the candidate, and reflected through 
the origin. This sample of candidates has been shown to be an unbiased estimator of the 
combinatorial background [1^, |1^ and to reliably model the shape of the background. The 
contribution from gluon splitting is estimated and subtracted from the sample based on the 
OPAL measurement of the multiplicity of such events in hadronic Z° decays [|18|. For the 
cuts used in this analysis, g — >■ cc events contribute (1.1 ± 0.4)% to the signal. After all 
corrections, and after combinatorial background subtraction, (15784 ±99) D*"*" mesons are used 
in the subsequent analysis. The error quoted is the statistical uncertainty of the combinatorial 
background subtraction. 
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Figure 1: Distributions of the mass difference AM = Mj)*+ — M-^o reconstructed in the four different 
D*+ channels. The arrows indicate the range in AM considered as signal. The background estimator 
distributions are superimposed, normalised to the signal distribution at large values of AM indicated 
by the cross-hatched area. Note that the significant tails in the AM distribution above the expected 
signal, particularly in (c) and (d), are caused by partially reconstructed D*"*" mesons, and is properly 
treated by the background estimator (see text). 
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Error source 


Variation 


Error 


Total uncorrelated 


±0.021 


B mixing 

Fragmentation modelling 
Gluon splitting 


Xeff : +11% 

(xs)b : +0.008, (xs)d : +0.009 
ng_cc : +21% 


-0.002 
+0.004 
-0.002 


Total 


± 0.022 



Table 2: List of the systematic errors on the charm fraction f^* in the D*+ sample. The top part 
of the table contains that part of the error which is uncorrelated with the systematic error associated 
to the reconstruction of leptons in the D*"*" sample. The signs given for the errors in the lower part 
indicate the direction in which the result changes for a change of the relevant variable by the amount 
and direction indicated in the middle column. 



The remaining two sources of D*+ production, —>■ cc ^ D*+X and ^ hh —>■ D*+X, are 
separated by applying three different flavour tagging methods, based on lifetime information, jet 
shapes and hemisphere charge information, as described in [r3|. Combining all D*+ channels, 
the overall charm fraction is determined to be: 

= 0.774 ± 0.008 ± 0.022, (7) 

where the first error is statistical and the second systematic. The dominant systematic errors 
are from the estimation of the background in the D*"^ sample, and from modelling the charm 
physics parameters used in the flavour separation. A breakdown of the systematic error into 
its components which are relevant for this analysis is given in table |^. The errors are split into 
two groups: one group which is uncorrelated to errors encountered when identifying leptons 
in this sample of events, and a second group of correlated errors. In the latter case the errors 
are signed indicating in which direction the result changes if the underlying physics variable 
is changed in the direction indicated in the table. More details of the procedure and of all 
systematic errors are given in [O]. 



5 The D*+£- Sample 



The D*"*'£~ sample is found by searching the hemisphere opposite the identified D*"*" meson for 
a lepton with a charge opposite to that of the D*"*" candidate. Electrons are identified using a 
neural network technique . The network used in this part of the analysis is slightly simplified 



compared to the one used in |]14|, using only 6 inputs, 8 nodes in one hidden layer, and one 
output. The input variables are 

• the difference between the measured specific energy loss, dE/dx, and that expected for 
an electron, divided by its expected uncertainty; 

• the experimental uncertainty on dE/dx; 

• E/p, the energy of the electromagnetic cluster associated with the track inside a cone 
with a half opening angle of 30 mrad, divided by the measured track momentum; 

• the number of electromagnetic blocks in the cluster; 
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Efficiencies for 


Source 


Electrons 


Muons 




0.302 ± 0.007 


0.433 ± 0.008 




0.305 ± 0.015 


0.305 ± 0.015 




0.222 ± 0.013 


0.346 ± 0.015 




0.213 ± 0.031 


0.315 ± 0.036 



Table 3: Efficiencies to reconstruct an electron or a muon opposite a D*^ meson separately for the 
different sources after applying all cuts. The errors quoted are purely statistical. 



• the momentum of the track; 

• the polar angle, | cos6'|, of the track. 

All variables are well modelled in the Monte Carlo simulation, thus ensuring a reliable calcula- 
tion of the selection efficiency. 

Muons are identified based on the of the matching between track segments in the central 



tracking chambers and in the muon chambers |]l5|. In addition the specific energy loss, dE'/dx, 



has to be compatible with that expected for a muon at the measured momentum. 

To reduce systematic uncertainties, electrons are reconstructed only in the central part of 
the OPAL detector, |cos^| < 0.715, while muons are required to satisfy |cos^| < 0.9. To 
increase the purity of the electron and muon samples, candidate tracks must have momenta 
greater than 2 GeV/ c. Events from bottom decays are suppressed by selecting only candidates 
with pt < 1-2 GeV/c for both electrons and muons, where the transverse momentum, pt, is 
measured with respect to the axis of the jet containing the lepton candidate, including the 
lepton candidate itself in the jet-axis calculation. After all cuts, a total of 661 electron and 
1045 muon candidates are selected. 

The efficiency to select a lepton is calculated using Monte Carlo simulations. It is calculated 
from events where a D*"*" meson is reconstructed in one hemisphere, and a lepton in the other, so 
that possible correlations between both hemispheres are taken into account. The e^^i efficiencies 
are found to be 

ec^e = 0.302 ± 0.007 and e^^^ = 0.433 ± 0.008, (8) 

where the errors are due to the finite Monte Carlo statistics. A list of all efficiencies, including 
those for leptons in bottom events, is given in table ^ The overall difference in the efficiencies 
for muons and electrons is mostly due to the larger range of cos 6' used for the muons. The 
ratio ey^^g/ec^i is larger for electrons than it is for muons, because the electron identification 
algorithm depends more strongly on pt than the muon identification does, the former being 
more efficient at large pt. Since h ^ £ events have on average a larger pt, electrons are found 
with larger efficiency in b — ^ events. The pt and p distributions of the selected candidates are 
shown in figures ^ and respectively. The distributions of the wrong sign candidates are 
superimposed. 

5.1 Combinatorial Background Estimation 

Background in the D*~^i~ events is estimated from the data with the help of the wrong sign 
D*+£+ sample. Subtracting the number of events found in the wrong sign sample from the 
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Figure 2: Transverse momentum spectrum (a) and momentum spectrum (b) of the selected lepton 
candidates. The arrow in (a) indicates the position of the pt cut. The hatched distribution is the 
background estimated using tlie wrong sign event sample. Composition of the pt spectrum in the 
Monte Carlo for right sign (c), and for wrong sign events (d). 
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Source 


Branching 

ratio ref. 


Sample co 
Electrons 


mposition 
Muons 


c-^ ^ 

b ^ c ^ £ 
b^c^£ 

b ^ J/^ ^ £+r 

Non-prompt 
Hadrons ^ 


0.1099 ± 0.0023 [g] 
0.0780 ± 0.0060 [|] 
0.0130 ± 0.0050 [IC] 
0.0045 ± 0.0007 |] 
0.0007 ± 0.0001 [|] 


0.604 ± 0.019 
0.083 ± 0.011 
0.135 ± 0.014 
0.030 ± 0.006 
0.003 ± 0.002 
0.003 ± 0.002 
0.063 ± 0.009 
0.079 ± 0.009 


0.580 ± 0.010 
0.057 ± 0.007 
0.098 ± 0.010 
0.017 ± 0.004 
0.005 ± 0.004 
0.001 ± 0.004 
0.089 ± 0.008 
0.153 ± 0.011 



Table 4: Semileptonic branching ratios as given in [T^ and [|l| and composition of the leptons in the 
D*"*" sample, as found in the Monte Carlo. Note that the sample composition is given for information 
purposes only, and is not used in the actual analysis. 



number of events found in the right sign sample gives an estimate of the number of c — £ 
decays, with only a small contribution remaining from background events (see equation 
The compositions of the right sign and of the wrong sign samples are shown in figures ^c and 

The subtraction of the combinatorial background relies on the assumption that these events 
are equally distributed between the right and the wrong sign sample, namely that = N^^. 
This subtraction procedure requires no explicit knowledge of the hadronic contamination in the 
lepton sample, since it is subtracted together with the wrong sign events. In figure |^, the 
shape of the pt distribution of the combinatorial background in the right sign sample, N^g^, is 
compared to the combinatorial component in the wrong sign sample, N^^. Good agreement 
is observed for the fraction of events below the applied cut of 1.2 GeV in the right and in the 
wrong sign combinatorial background. The shapes are slightly different which is attributed to 
different contributions from bottom events to both samples. However since only the overall 
number of events is needed in this analysis the difference has a very small influence on the final 
result. 

Monte Carlo studies show that the assumption = is not entirely correct, since a 
particular class of events, accounting for less than 10% of the background, is found more often 
in the right sign sample than in the wrong sign sample. These events consist of a partially 
reconstructed D*"*" meson opposite a correctly identified lepton with the correct charge correla- 
tion. The number of such events found in the wrong sign sample amounts to only 55% of the 
number of the same type of events found in the right sign sample. The total number of these 
events in the right sign sample has been measured in |[1^ from data. Relative to the combina- 
torial background, they account for (8.5 ±2.2)% of the total right sign sample. The background 
subtracted sample is therefore corrected for the fraction of these events, namely by +45% of 
the (8.5 ± 2.2)%. This corresponds to a background charge asymmetry of AA^bgd = +(43 ± 11) 
events, where the error is dominated by the fraction of such events measured in the data. An 
additional modelling error of 50% of this correction is applied, as will be discussed in section |^. 
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Figure 3: Comparison of the right sign combinatorial background (points with error bars) and the 
wrong sign combinatorial background component (line histogram) in the simulation. 



5.2 Estimation of the Bottom Background 

The sample of tagged D*^i~ events has a charm purity of about 60%. Non-leptonic background 
accounts for 8% of the electron candidates and 15% of the muon candidates. The rest consists 
of correctly identified leptons from a number of different sources. The sample composition as 
determined from the Monte Carlo simulations is shown in table H, and illustrated in figures 
^ and l^d. Since the charges of the D* and the i should be opposite in the D*+^~ sample, 
any effect which influences the charge correlation between the two hemispheres influences the 
flavour composition. The most important of these is B-meson mixing. If mixing has occurred 
in either hemisphere, the charge correlation between the primary quark and the corresponding 
D*"*" meson is changed. The total probability in b-events that mixing has changed the charge 
correlation is given by 



XD*+e- 



Xd*+(1 - Xe) + X^l - Xd*+) 



(9) 



where Xd*+ ? Xe are the effective mixing parameters applicable to the D*+ and the lepton, re- 
spectively. These effective mixing parameters depend on the fractions of and B^ mesons in 
the sample under consideration, and on the mixing in the B[J and the B^ system. The average 
mixing in the B° system is measured to be Xd = 0.175 ± 0.016 The LEP combined lower 
limit for Bg mixing given in corresponds to a lower limit on Xs of 0.49 at 95% confidence 
level. In this analysis, Xs is varied between 0.49 and the maximum value of 0.50. 

Most D*"*" mesons in bb events originate from decays of the B[| meson. In |19| this 

fraction has been determined to be rj*^ = O.Sllgji- The fraction of D*+ that come from 
Bg mesons has been estimated to be r^*^ = 0.043 ± 0.039 |]TU[. The effective mixing in the 
hemisphere containing the D*"^ meson is therefore 



Xd* 



■Xd + r, 



■ Xs 



0.163: 



-0.025 
-0.030 



(10) 



The fraction of leptons produced in decays of B|j and B^ mesons is determined from the fractions 
of weakly decaying B-hadrons in Z° — > bb events by weighting with the lifetimes of the B- 
hadrons species [ffl]. This is done in order to correct for the different semileptonic branching 
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ratios and leads to the values = 0.399±0.023 and rf = 0.1181q;o2o, respectively. The effective 
mixing parameter is 

Xi- = r^ Xd + ri ■ Xs = 0.128 ± 0.012 . (11) 

In addition, D*"*" mesons with the wrong sign can be produced in bottom decays, where a c 
quark is produced in the decay of a virtual W. This can be expressed in terms of a mixing-like 
parameter As in [^, a value of Cd = 0.025 ±0.025 is used. The effective mixing parameter 
is then 

Xeflf = XD'+e-i'^ - Cd) +Cd(1 - Xi)*+e-) , (12) 

neglecting terms which are quadratic in either XT>*+i- or Cd- The effective mixing parameter 
for the D*+£~ sample is finally estimated to be 

XeS = 0.2611^34- (13) 

In total, the contribution to the background from bottom events is calculated according to 
equation ^, using the efficiencies listed in table |^, the branching ratios given in table ^ and 
the effective mixing parameter determined above. The total contribution amounts to AA^b = 
— (53 ± 7) events. 



6 Results 

The number of D*^i~ combinations in charm events is determined according to equation |^. 
The background subtracted momentum and transverse momentum spectra for electrons and 
muons are shown separately in figure ^. The distributions are further corrected for the effects 
of mixing and for the charge asymmetry in the background, as described in section ETT. 



The total number of leptons from charm hadron decays is iVp,+ g- = 378 ±31 and A^,+ ^- = 
476 ± 40, respectively. The quoted error is purely statistical. In table ^, a summary of the 
selected events in each sample is shown. Combining these measurements with the total number 
of selected D*"*" mesons, Nd*+ = 15784 ± 99, the appropriate charm fraction, and the lepton 
efficiencies, the inclusive semileptonic branching ratios of charm hadrons in cc events are 

determined to be 

B(c ^ e) = 0.103 ±0.009 and B(c ^ /x) = 0.090 ± 0.007 . 

Here, the quoted errors are only statistical. The semileptonic branching fraction of charm 
hadrons derived from these individual results is 

B(c^i) = 0.095 ±0.006. 



7 Systematic Errors 

In this section, the different sources of systematic errors are discussed. A breakdown of all 
errors considered is summarised in table ^. All errors in this section are given relative to the 
inclusive B{c ^ i) branching ratio, if not otherwise stated. 
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Figure 4: Momentum spectra after background subtraction for electrons (a) and muons (b), and 
transverse momentum spectra for electrons (c) and muons (d). Points are data, the line histogram is 
the Monte Carlo prediction. Both data and Monte Carlo include the residual background contributions 
from bottom events and from the background charge asymmetry. 
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Sample 


D*e 


D> 


Right- sign 
Wrong-sign 


661 
305 


1045 
558 




378 ± 31 


476 ± 40 



Table 5: Summary of selected events in each sample, and number of events after background 
subtraction. 



Modelling errors 

— c — s> £ modelling: The momentum spectrum of leptons in c — ^ ^ decays is described 



by the ACCMM model. Using the range of parameters recommended in |jTO| this 
corresponds to an error of ^3 3% of the charm semileptonic branching ratio. The 
size of the error is largely dependent on the momentum cut used in the identification 
of electrons and muons. 

h ^ i modelling: The momentum distribution of the leptons in bottom decays influ- 
ences the tagging efficiencies. Following the recommendations in this has been 



studied by reweighting the lepton spectrum in the Monte Carlo simulation to differ- 
ent theoretical models, with ranges of parameters chosen such that the experimental 
errors are covered. The ACCMM model is used to obtain the central value, and 
the ISGW and the ISGW** p2[ models are used for the iblcr variation around the 
central value, and the efficiencies are recalculated. The errors found are 0.1%. 

h —>■ c i modelling: The ACCMM model is used to describe the momentum 



spectrum of cascade h c ^ i decays, as suggested in |T0|. Three different sets of 

parameters are proposed to cover the experimental uncertainties in the momentum 

spectrum. The lepton efficiencies are recalculated. The errors for the final result are 

+0.20/ 

-0.1 /o- 

Fragmentation modelling: The fragmentation parameters in the Monte Carlo simu- 
lation have been varied to change the mean scaled energy of weakly decaying bottom 
and charm hadrons around their experimental values of {xe)b = 0.702 ± 0.008 and 
{xe)b = 0.484 ±0.009 respectively |T0|. This study is done using the Peterson frag- 
mentation model 1^. This results in an error of ±0.6%. In addition, the Peterson 
fragmentation model has been replaced by the Collins and Spiller fragmentation 
model p3[ and by the Kartvelishvili fragmentation model |2^. The parameters for 



these models have been adjusted to the same mean scaled energy as for the Peterson 
function. The largest deviation between the different models is used as a systematic 
error. Combined with the error using different parameters for the Peterson model, 
a total error of 0.9% is determined. 



B-physics 



B-meson mixing: The uncertainty due to mixing in the neutral B sector has been 
studied by varying the effective mixing parameter (see equation |T3[) Xcfr within its 
errors. An error of 0.8% is found. 

Branching ratios: The dependence on the branching ratios b — >■ £ and b ^ c — ^ has 
been investigated by varying them within their experimental errors. Mean values 
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and errors used are given in table ^ An error of 0.8% is found. A breakdown of the 
error into the different channels contributing is given in table |^. 

• Particle identification 

— Electron identification: The efficiency to identify electrons is calculated in the Monte 
Carlo. The two variables which mainly determine the performance of the neural 
network are the specific energy loss, dE/dx, together with its error, and the ratio 
E/p. Both variables are compared between Monte Carlo and data using different 
samples of identified particles. The dE/dx measurements are calibrated in data 
using samples of inclusive pions at low momenta and electrons from Bhabha events 
at 45 GeV/c. The quality of the calibration is checked with a number of control 
samples, mostly pions from Ks decays and electrons from photon conversions. The 
deviation between the mean dE / dx measured for these samples in the data, and the 
mean dE/dx in the Monte Carlo, is below 5%. Similarly, the resolution of dE/dx 
is studied in these samples, and is found to be described in the Monte Carlo to 
better than 8%. The total error from these two effects is found by varying both 
simultaneously, and is ±2.5% Note that for this analysis, no explicit knowledge of 
the hadronic background in the sample of lepton candidates is needed, since it is 
subtracted using the wrong sign sample. 

A similar study has been performed for the next most significant input variable, 
E/p. The E/p resolution in the Monte Carlo is about 10% worse than in the data. 
The Monte Carlo has been reweighted to the data, and the full difference is used as 
an estimate of the error, resulting in a variation of the efficiency of ±2.7%. 
No significant contributions to the error are found from the other input variables of 
the network. The error related to them is estimated from the statistical precision 
of these tests, which is less than 1% of the efficiency. In total, an error of ±4% is 
assigned to this source. 

— Muon identification: The systematic error of the muon identification efficiency is 
evaluated using a method similar to that described in [0. The muon detection 
efficiency is compared between data and Monte Carlo using various control samples, 
namely Z° ^^fi~ events, and muons reconstructed in jets. Without using dE/dx 
information, an error of ±2% is found. The influence of the dE/dx selection cut 
on the muon ID is studied in the same way as described for electrons. The mean 
dE/dx for muons in Z° —>■ fi'^fj^~ events is observed to be shifted by approximately 
15% of the resolution in dE/dx with respect to the theoretically expected value. A 
very similar shift is observed in the Monte Carlo, both for muon pairs and for muons 
identified inclusively in jets. An error of 5% is used. 

The dE/dx resolution is studied in the data, and is found to be modeled by the 
Monte Carlo to better than 5%. The final error assigned to the efficiency of muon 
identification is ±3.0%. 



• Internal sources 



Flavour separation: The errors of the flavour composition on the D* 
timated in 
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sample es- 

are used to calculate the corresponding error of the semileptonic 
branching fractions. A breakdown of the total error into sources correlated and un- 
correlated with the reconstruction of leptons in the D*"*" sample is given in table |^, 
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and is taken into account in calculating the final error. The uncorrelated error is 
2.8%; the total correlated error is 0.9%. 

Background charge asymmetry: The correction applied to the background-subtracted 
sample of D*~^i~ events is calculated based on the measured fraction of events con- 
tributing, and on the charge asymmetry, which comes from Monte Carlo simulation. 
For the former, the statistical error of the measurement is used as a systematic uncer- 
tainty, translating into an error of 1.1%. The Monte Carlo prediction of the charge 
asymmetry is conservatively varied by ±50% of its value. The final error from this 
is 2.5%. 

Background estimation: The background in the D*+£~ sample is estimated from the 
wrong sign sample. The number of combinatorial background events, corrected for 
mixing and for the effects of the background charge asymmetry, is compared with the 
expected number of combinatorial background events in the Monte Carlo simulation. 
Within the statistical precision of this test good agreement is found. The statistical 
error of this test is used as a systematic uncertainty, resulting in an error of 1.8%. 
The influence of the pt cut on the background is studied by comparing the shapes of 
the background between data and Monte Carlo. The data spectrum is reweighted 
to the Monte Carlo one, and the number of background events is recalculated. The 
resulting difference is used as a systematic error of 0.3%. The final error assigned is 
1.8%. 

Detector modelling: The influence of the detector resolution on the tagging vari- 
ables is studied in Monte Carlo simulations by varying the resolutions in the central 
tracking detectors by ±10% relative to the values that optimally describe the data. 
The analysis is repeated and the efficiencies are recalculated. The error is 1.1%. 
The calculation of the efficiencies relies on the correct modelling of the detector 
acceptances, in particular in cos^. This has been tested by reweighting the cos 6 
distribution of D*"*" candidates as found in the Monte Carlo simulation to that re- 
constructed from data, and repeating the analysis. This changes the result by 0.3%, 
which is used as a systematic error. The total the error due to detector modelling is 
1.2%. 

Gluon splitting: Gluon splitting into a pair of heavy quarks can produce D*"*" mesons 
which might contribute to the sample of selected events. This contribution is found 



to be (1.1 ± 0.4)%. It is based on the OPAL measurement of gluon splitting [|TS 
and Monte Carlo simulation to determine the selection efficiency. The total number 
of D*"*" mesons is corrected for this effect. The uncertainty of this number is used 
as a systematic error of 0.4%. Similarly, leptons can be produced in gluon splitting 
events. The contribution to the sample is found to be (0.2 ± 0.1)%, which results 
in an error of 0.1%. According to these studies the total systematic uncertainty is 
0.4%. 

• Monte Carlo statistics 

— Monte Carlo statistics: The efficiencies to identify a lepton in the D*"*" sample are 
calculated from the Monte Carlo with limited statistical precision. The error from 
this source amounts to 1.5%. 

A complete list of systematic errors is presented in table |^ for B(c e), B(c fi), and 
B(c — > i). Except for the error from Monte Carlo statistics and the lepton identification errors. 
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Soil TOP 


B(c ^ e) 


B(c ^ ii) 


B(c ^ i) 


Modellmg 




c — > £ model 


+0.0057 
-0.0034 


+0.0050 
-0.0030 


+0.0053 
-0.0031 


b — > £ model 


0.0001 


0.0001 


0.0001 


L) C lliULlt:! 


0.0002 


0.0002 


0.0002 


Fragmentation modelling 


0.0010 


0.0008 


0.0009 


Total modelling 


+0.0058 
-0.0035 


+0.0051 
-0.0032 


+0.0054 
-0.0032 


B physics 




B-meson mixing 


0.0009 


0.0008 


0.0008 


B(b £) 


0.0003 


0.0002 


0.0003 


B(b -> c £) 


0.0007 


0.0006 


0.0006 


B(b ^c^i) 


0.0005 


0.0004 


0.0005 


Hyo — > T — > I) 


0.0001 


0.0001 


0.0001 


Bfb ^ J/* ^ £'^£~) 


< 0.0001 


< 0.0001 


< 0.0001 


Total B physics 


0.0013 


0.0011 


0.0012 


Particle ID 




Electron identification 


0.0041 




0.0017 


Muon identification 




0.0027 


0.0015 


Internal sources 




Flavour separation (uncorr.) 


0.0028 


0.0024 


0.0026 


Background charge asymmetry 


0.0026 


0.0023 


0.0024 


Background estimator 


0.0019 


0.0016 


0.0017 


Detector modelling 


0.0012 


0.0011 


0.0011 


Gluon splitting 


0.0005 


0.0004 


0.0004 


Monte Carlo statistics 


0.0024 


0.0016 


0.0015 


Total error 


+0.0088 
-0.0075 


+0.0072 
-0.0060 


+0.0074 
-0.0060 



Table 6: List of systematic errors contributing to B(c ^ e), B(c /x) and B(c — >■ ^). A detailed 
explanation of the different errors can be found in the text. 
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all errors from a given source are assumed to be fully correlated between the electron and the 
muon results. 

To check the stability of the results, the analysis is repeated with different selection cuts 
for the leptons. Consistent results are found if the momentum cut is raised from 2 GeV/c to 
3 GeV/ c both for electrons and muons, if the transverse momentum cut is removed, if the muon 
selection is repeated using muons in the central part of the detector only, and if the muon 
selection is done without using the dE/dx selection cut. 



8 Conclusions 

A measurement of the inclusive charm hadron semileptonic branching fractions in Z° — > cc 
events, B(c e) and B(c — /i), has been presented. The identification of Z° — >■ cc events is 
based on the reconstruction of D*"*" mesons. The semileptonic branching ratios are measured 
by reconstructing leptons in the charm-tagged sample and are found to be 

B(c ^ e) = 0.103 ± 0.009l!^;[![J^ and B(c ^ /i) = 0.090 ± 0.007l|];E^|]^ , 

where the first error is in each case statistical and the second systematic. Combining the two 
measurements while taking correlations into account, gives 

B(c ^i) = 0.095 ± 0.006^°:°°^ . 

This result agrees well and is competitive with the most recent published measurement at lower 
energies of B(c ^ £) = 0.095 ± 0.009 ||. 
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